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Abstract. V889 Herculis is one of the brightest single early-G type stars, a young Sun, that 
is rotating fast enough (P ro t = 1.337 days) for mapping its surface by Doppler Imaging. The 
10 FOCES spectra collected between 13-16 Aug 2006 at Calar Alto Observatory allowed us to 
reconstruct one single Doppler image for two mapping lines. The Fe 1-6411 and Ca 1-6439 maps, 
in a good agreement, revealed an asymmetric polar cap and several weaker features at lower 
latitudes. Applying the sheared-image method with our Doppler reconstruction we perform an 
investigation to detect surface differential rotation (DR). The resulting DR parameter, 50. /Q, ~ 
0.009 of solar type, is compared to previous studies which reported either much stronger shear 
or comparably weak DR, or just preferred rigid rotation. Theoretical aspects are also considered 
and discussed. 
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1. Introduction 

The primary source of all the manifestations of the solar activity is the magnetic dy- 
namo. However, the Sun is the only star, for which one of the conditions of a working 
dynamo, the surface differential rotation can be measured directly, so far. Today's ad- 
vanced techniques can measure the solar surface and radial rotation profiles with high 
accuracy, but at the dawn of extensive solar observations, tracing of sunspots at different 
latitudes was the only tool to measure the surface shear on the Sun (e.g. Maunder & 
Maunder 1905). Similarly, on stars the detection of surface differential rotation is still 
a challenging observational task. Direct tracing of starspots is the usual way to observe 
surface shear, however, it requires a reliable surface reconstruction technique such as 
Doppler imaging. 

Detecting surface differential rotation (DR) on stars of different types can provide es- 
sential observational input for the theoretical understanding of DR and dynamo. Indeed, 
investigating solar analogues, such as our target, V889Her (HD 171488), provide us with 
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Table 1. Physical parameters of V889Her taken from Frasca et al. (2010) 



Paramete 



Vali 



Spectral type 
T cf f (K) 
P„t (days) _ 
v sin i (km s — x ) 
Inclination (°) 

Microturbulcncc 



(kn 



Macroturbulcnce £r — £t 



G2V 

5750±130 

1.33697±0.0002 

37.1±1.0 

60±10 

4.30 ± 0.15 

1.6 

3.0 



useful information for further development of solar dynamo theory, which is still far from 
a thorough understanding. 

The fast-rotating solar-type V889Her is probably the brightest single early-G type star 
that is rotating fast enough for mapping its photosphere by means of Doppler imaging. 
The star is a member of the Local Association, a stream of young stars with ages rang- 
ing from 20 to 150 Myr (Montes et al. 2001). Strassmeier et al. (2003) made the first 
photometric and spectroscopic study of V889 Her and presented the first Doppler recon- 
struction showing a large polar spot and additional high-latitude features. Subsequent 
works by Marsden et al. (2006) and Jeffers & Donati (2008) based on Zeeman-Doppler 
imaging have also shown polar spottedness and found a strong solar-type differential 
rotation. However, Jarvinen et al. (2008) reported instead a much weaker surface shear, 
while Huber et al. (2009) preferred rigid rotation. In a comprehensive study of the chro- 
mospheric and photospheric activity, including Doppler imaging, Frasca et al. (2010) 
concluded that V889 Her is a G2V type star. Hereafter, for further investigation of DR 
we adopt the reviewed astrophysical parameters given there (see Table Q]) . 



2. Sheared-image reconstruction 

In order to detect differential rotation, the Doppler reconstruction of Frasca et al. 
(2010) is repeated, but instead of assuming rigid rotation, latitude-dependency is al- 
lowed in the rotation profile. This technique known as 'sheared-image method' (Donati 
et al. 2000) incorporates a predefined solar-type differential rotation law in the Doppler 
imaging process. In practice, the image shear expressed by a = 5Q/Q is kept fixed, 
and individual Doppler reconstructions are done over a reasonable range of values in 
the P ro t — ol parameter plane. Each of the resulting Doppler maps has a % 2 -value which 
describes the goodness-of-fit to the observed spectral line profiles. The best combination 
of Prot and a is then marked by the least x 2 value. The ability of this method was 
demonstrated also on artificial data by Kovari & Weber (2004). 



Fel-6411 V889Her 2006.7 _ Cal-6439 V889Her 2006.7 




Longitude (deg) Longitude (deg) 



Figure 1. Doppler maps for the Fe 1-6411 lines (right) and for the Ca 1-6439 lines (left), taken 
from Frasca et al. (2010). The arrows below indicate the phases of spectroscopic observations. 
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Figure 2. Resulting \ maps of the sheared-image method over the P ro t — a parameter plane. 
The best fit is found at a = 0.003 ± 0.004 and P ro t = 1.3395 ± 0.0027 days for the Fe 1-6411 
reconstructions {top left panel) and at a = 0.010 ±0.002 and P rot = 1.3350 ± 0.0052 days for the 
Ca 1-6439 reconstructions {top right panel). Bottom panel shows the averaged \ 2 map, where 
minimum is found at a = 0.009 ± 0.006 and at P rot = 1.3395 ± 0.0028 days. 



For the Doppler reconstructions we used the Doppler imaging code TempMap by Rice 
et al. (1989). Doppler imaging was performed for the Fei-6411, and Cai-6439 mapping 
lines. The individual reconstructions shown in Fig. [1] revealed similar spot distributions, 
i.e., mainly cool polar spots with temperature contrasts of up to «1500K with respect to 
the unspotted surface of «5800K. Some low- latitude features are also recovered, however, 
with significantly weaker contrast ranging from «300K (Ca 1-6439) to a maximum of 
«500K (Fei-6411). 

Fig. [2] displays the x 2 -maps f the trial-and-error process over the P r ot — ct parameter 
plane. The best fit minima have similar locations for the two independent image recon- 
structions. We find £1(9) = 4.69 — 0.014 sin 2 9 rad/day differential rotation law for the 
iron line and il{9) = 4.71 — 0.047 sin 2 9 rad/day for the calcium line, while the combined 
X 2 map suggests £7 oq = 4.69 rad/day and Si} — 0.042 rad/day. 



3. Discussion 

Using the averaged x 2 landscape of the two independent line reconstructions in Fig. m 
we find a = 0.009 or equivalently SVl = fi oq — n po \c = 0.042 rad/day. Accordingly, the 
time the equator needs to lap the pole by one full rotation is about 150 days, i.e., of the 
same order of the solar value. Observations and theoretical model calculations indicate 
that differential rotation is firmly influenced by stellar temperature (cf. e.g., Barnes et al. 
2005, Reiners 2006, Kitchatinov & Riidiger 1995). A recent numerical model developed by 
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Kitchatinov & Olemskoy (2010) predicts SSI w 0.075 rad/days for a T off = 5800 K dwarf, 
which yields a « 0.016 at the angular velocity of V889Her. This shear is much weaker 
than SSI « 0.4 — 0.5 rad/day obtained by Marsden et al. (2006) and by Jeffers & Donati 
(2008), respectively, both from Zeeman Doppler imaging. Those values are among the 
largest ones measured only for very few stars of mainly F-type (Reiners 2006). Such a large 
value does not fit the power law of Barnes et al. (2005) which suggests SSI « 0.12 rad/day 
for a G2-type star. Moreover, Jarvinen et al. (2008) argues for a substantially weaker 
differential rotation. Indeed, Huber et al. (2009) preferred solid rotation, though, due to 
the rather low data quality, their results does not exclude a weak surface shear such as 
that was suggested by Jarvinen et al. (2008) and by us. 

All in all, the differential rotation derived from our sheared-image process for the G2- 
star V889Her is below but close to the one from the empirical relationship in Barnes 
et al. (2005) and seems to support the theoretical model predictions (Kiiker & Riidiger 
2005, Kitchatinov & Olemskoy 2010) as well. On the other hand, such a weak differential 
rotation is inconsistent with the high surface shear derived from Zeeman Doppler imaging 
(Marsden et al. 2006, Jeffers & Donati 2008). This difference might partially be explained 
by temporal variations in differential rotation, however, in the case of V889Her Jeffers 
et al. (2010) found no evidence for such a temporal change. If so, other effects such as 
rapid spot evolution should also be considered. In addition, the very different methods 
and data applied to detect the surface shear could have different, sometimes depreciated 
or unrecognized error sources, which can also explain such a disagreement (see also T.A. 
Carroll's paper in this volume). For instance, when various spectral features used for 
imaging correspond to different surface layers, the detected shear would diverge. 
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